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abstract: A time-discretization that preserves the super-integrability of the Calogero model is 
obtained by application of the integrable time- discretization of the harmonic oscillator to the 
CS| ■ projection method for the Calogero model with continuous time. In particular, the difference 
equations of motion, which provide an explicit scheme for time-integration, are explicitly presented 
for the two-body case. Numerical results exhibit that the scheme conserves all the(= 3) conserved 
quantities of the (two-body) Calogero model with a precision of the machine epsilon times the 
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I Introduction 



Numerical analysis of dynamical systems has great importance and a wide variety of applica- 
tions in science and engineering. The elaboration of schemes for numerical analysis has a long 
and continuous history of studies as well as a rich accumulation of techniques. Whenever one 
applies numerical analysis to the equations of motion, one must discretize the time-evolution of 
the dynamical system that is originally described by differential equations because of the lack 
of the notion of infinity in numerical analysis. This leads to difference equations, which do not 
usually describe the same dynamical system as the original one: time-discretization is typically 
accompanied by modification of the original system, which may cause a significant difference in 
the behavior of the solution from that of the original system, particularly after integration over 
a long period. Understanding and controlling such modifications are thus important in the quest 
for more accurate long-time integration in numerical analysis. 

The symplectic integration method (see, for instance, refs. [U 121 El), or the symplectic inte- 
grator, is one of the time-discretizations that was invented in such a quest. Given a Hamiltonian 
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, it is designed so that its one-step(= r) time evolution gives the exact one-step time evolution 
of a modified Hamiltonian H := H^ ' + tH^ + t 2 H^ + • • • . Since the modified Hamiltonian 
is conserved by the flow of the symplectic integrator, the fluctuation of the value of the original 
Hamiltonian H^°\ i.e., the total energy of the system with continuous time, is bounded, which 
is far more favorable than unbounded increase or decrease of the total energy that one usually 
observes in other non-symplectic discretizations. This is the reason why the symplectic integration 
method shows better accuracy even after long-time integration. 

However, even the symplectic integration method does not usually have modified constants of 
motion for all the constants of motion of a system, which often causes secular increase or decrease 
of the values of the constants of motion after long-time integration. For example, non-existence 
of the modified constants of motion is proved for the two super-integrable models discretized by 
the symplectic integrator: the two dimensional Harmonic oscillator with integer frequency ratio 
(except for the isotropic case) and the two-dimensional Kepler problem [HI5]. The orbits generated 
by the symplectic integrator are not closed, though those of the exact analytic solutions are closed, 
indeed. 

Thus the following question naturally arises: are there any discretization schemes that preserve 
the (super-)integrability of (super-)integrable models? Actually, an extensive collection of the 
known integrable discretizations of integrable models is now available in the monograph [Oj that 
was published in recent years. However, the super-integrable discretization of super integrable 
models has not been studied that far yet. Quite recently, an affirmative answer to the question on 
the super-integrable discretization is shown for the Kepler problem for two and three dimensional 
cases IS], where the integrable discretization of the harmonic oscillator |U| plays an essential 
role. The above super-integrable discretization conserves all the constants of motion, i.e., the 
Hamiltonian, the angular momentum and the Runge-Lenz vector, and generates a sequence of 
discrete points on the orbit of the exact analytic solution of the Kepler problem. And, of course, 
the orbit with the eccentricity less than unity is closed. 

The purpose of this paper is to present a super-integrable discretization of the Calogero 
model HQlin], 

i=l i,j=l V 1 31 

The real- valued quantities pi, Xi, u and a in the Hamiltonian are the canonical momentum and co- 
ordinate of the i-ih particle, the strength of the external harmonic confinement and the interaction 
parameter, respectively. The Calogero model (u ^ 0) with N degrees of freedom (corresponding 
to the iV-body case) is maximally super-integrable in the sense that it has 2N — 1 constants of 
motion which are independent of each other ^2]- The super-integrable structure of the Calogero 
model is built up by the Lax formulation fSlE], with which the eigenvalue problem of an oscillat- 
ing Hermitian matrix |T3] is intrinsically involved. We shall present a discretization that preserves 
the above super-integrable structure of the Calogero model for the general iV-body case. It gives, 
in particular, the explicit form of the difference equations of motion of the Calogero model for the 
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two-body case that conserves all the three constants of motion as 



j 2 At 2 



Pl,n ~ —X 1<n At 



2 At 2 



1 



2a 



J 2 At 2 



-M Un At 



]_ , uj 2 At 2 Xl > n + -| 



+ 



A + p ltn = - 



1 



4 

j 2 At 2 



o 2 At 2 
4 



Ll + 



j 2 Ai 2 



CO 



1 + 
1 



u 2 A*2 



(x 1)Tl + £ 2 ,n) + 



1 



1 + 



M r;n At 

5^-(pi,« + P2,n)Atl + Ml n At 2 )At, 



+ 



1 



w 2 At 



2 



2a 



-M, n 



1.2) 



1 + 



^ 2 At 2 
4 



+ 



4 

2 Ai 2 



1 + 



j 2 At 2 



x hn At 



1 + 



j 2 At 2 



with 



A_|_X2,n — A + Xi ]n 



l<->2 



v/4a 2 At 2 + F n 2 
cu 2 At 2 i 



(Pl.n + £>2,n) - 1 
A+p 2 ,n = &+Pl,n 

M r , n 



M r , n At 

+ z 2 ,„)At| (Mi 2 n + M 2 n At 2 ) At, 



2a 2 



4 



4a 2 At 2 + F n 2 + F nV /4a 2 At 2 + F n 2 ' 

(aCl,n - a; 2,n) 2 + (Pl,n ~ P2,n)(^l,n ~ &2,n)At, 



F ■ = 
where a; i)n , 

discrete time. The symbol A + denotes the advanced time-difference defined by 
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for an arbitrary variable A n (e.g., A n 

The paper is organized as follows. In section m we present a brief summary of the projection 
method ^H] (for review, see refs. [HilEllIE!, for example), which gives a solution to the initial 
value problem of the Calogero model with continuous time. The Lax equations for the Calogero 
model provide a map from the Calogero model into the matrix-valued harmonic oscillator and 
also play an essential role in our discretization. Applying the integrable discretization for the 
harmonic oscillator 0, we discretize the projection method in section ITTT1 A discrete analogue of 
the Lax equations, which we call the dLax equations, is derived as a natural consequence of the 
discretization. In section II VI we present the explicit forms of the dLax equations of the Calogero 
model for the two-body case, which are equivalent to the difference equations of motion (jl.2j) . They 
provide an explicit scheme for the time-integration of the model. Numerical results obtained by 
our integrable discretization as well as by two other discretization schemes, namely, the symplectic 
Euler and energy conservation methods, are also presented. Sectionals dedicated to the summary 
and concluding remarks. 



II Projection Method 

In terms of the Lax pair for the Calogero-Moser model [T3] . 



L y (t) := Pi(t)<J fi + ^ -a; .(t) ^ ~ 

: = E , , t) 1" (t)f - , . t . 1° f a - u w = i.2,...,w, (2 ' 1, 

as well as a diagonal matrix 

:= diag(xi(t),x 2 (t), . . .,x N (t)), 
the canonical equation of motion for the Calogero model (jl. 1|) can be cast into the Lax form, 

^ = [L, M] - ^D, 

where [A, S] := — or equivalently, 
dL± 



[L , M] ± io»L , (2.2) 



df 
where 

L^t) :=L(t)±ia;D(t). (2.3) 

The Lax equations (j2.2j) allow the following relation between the products of L + and L~ and 
the Hamiltonian 

j t {{L + ) h (L-) mi • • • ) = \{L + ) h (L-r> ■ ■ • , M] + ico (h + ■ ■ ■ - (m, + • • • )) (L+)' 1 (IT)"* • • • , 
for any nonnegative integers li,mi, .... Iterated use of the above formula yields 

^(j[Tr(Y[(L + ) l ^(L-) mi ^ =^]^Tr(]^(L + )^(L-) m ^)Tr[]^(L + )' fc ^L-) mfc .SM 



i ij^k j 



k,l 

where 
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Thus the constants of motion of the Calogero model can be constructed by taking the trace of 
any products of L + and L~ of the following form, 

nTr(n(/- ) ; (/- )'"-)• if » ; »„• (2.4) 

i j i,j i,j 

which includes the constants of motion given in ref. jT^j. Considering the case a = where the 
matrices L^ 1 become diagonal, one can confirm the quantities given eq. (|2.4|) above include at least 
(and also at most) 2N — 1 constants of motion that are independent of each other. For example, 
one confirms 

N N 

C i: =Tr^TrL-=(^p„) 2 + - 2 (E^) 2 ? (2 

i=i i=i ^■ t> > 

h := TiL + L- = 2H, I 2 = Tr(L + ) 2 Tr(L~) 2 , 

are constants of motion of the Calogero model (jl.lj) that are independent of each other. Note that 
the relations satisfied by the Lax pairs (T^J 120] 

[L + , L-] = -2iu[L, D] = -2ua(T - E), T tj := 1, E i5 := 

M*[l, . . . , 1] = 0, [1, . . - , 1]M = 0, [^TM = MT = 0] U0) 

particularly the relations in the second line of eq. (|2.6|) which we call the sum-to-zero property, 
have played a crucial role in the quantum analogue of the Lax formulation. 

The initial value problem of the Calogero model can be solved by the projection method [To] , 
in which the Lax formulation presented above plays a crucial role. This method can be formulated 
in an analogous way to the Dirac picture in the time- dependent perturbation theory of quantum 
mechanics. Let us introduce L^, the L ± matrix in the "Dirac picture" by 

L±(f) =: ef^Lftf), L±(0) = L±(0). (2.7) 

Then the Lax equation (J2.2|) is rewritten as 

which has the same form as the Heisenberg equation in the Dirac picture with the time-dependent 
perturbation M(t) := M[Xi(t), ^(t), . . . , Xjy(t)) . 
The above equation allows the formal solution as 

L±(t) = U\0, ^(OMO, t), (2.8) 

where the time-evolution unitary matrix U(t', t) is given by the Dyson series of M(t), 

U(t', t) := J2 / d ^ / dt k-i ■ ■ dtxM^) ■ • ■ M(t k ^)M(t k ), (2.9) 



which has the semigroup property: 

U{0,t) := U[x(0),p(0);a,u)]t] 

= U[x(0), p(0); a, cj; ^[/[^(t'), p(t'); a,u;t- t'] = U(0, t')U(t', t), (2.10) 
C/t(0,t) = C/(t,0). 

Note that U(t',t) has a constant eigenvector *[!,..., 1] whose eigenvalue is unity, 



U(t', t)% . . . , 1] = *[1, . . . , 1], [1, . . . , l]U(t', t) = [1, . . . , 1] 
U(t',t)T = TU(t',t) = T. 



Xii) 



This property of ?7(t', t) is a consequence of the sum-to- zero property of the matrix M in eq. (|2.6J) . 

Substitution of the formal solution ()2.8|) into eqs. (|2.3j) and (j2.7j) gives the following solution 
of the initial value problem of the Lax equation (|2.2jl : 

JW = 1 /*(0, t ) ^^^ (0,«), 

(2.12) 

L( t ) = C /t(0, t ) e "' L+( ° )+ 2 e -' L - ( % (0, t ). 

e iaJ *^+(0) — e _iw *L _ (0) 

which means that the eigenvalues of the time-dependent Hermitian matrix — 

give the solution of the initial value problem of the Calogero model. The time-evolution unitary 
matrix U(0,t), which has been formally introduced as the Dyson series of M(t), is given here as 
the diagonalizing matrix. 

The unitary matrix £7(0, t) provides a map of the Calogero model into the matrix-valued 
harmonic oscillator. Let us introduce the matrices 



X ± (t) := U(0, t) = P(t) ± itoQ(t), 

P(t) ■= U(0, t)L(t)rf(0, t), Q(t) := U(0, t)D(t)U*(0, t). 

Substituting into the Lax equation (J2.2)) . we have 

dX ± 



[2.13) 



-- ±iuX ± } (2.14) 
which are equivalent to the equations of motion of the harmonic oscillator 

Thus we confirm that the Lax equations of the Calogero model ()2.2|) can be mapped to the 
equations of motion of the matrix- valued harmonic oscillator, eqs. ()2.14|) and ()2.15p . Though 
the Hermitian matrices P(t) and Q(t) can possess 2N 2 parameters for their initial values, their 
definitions (J2.1H)) introduce the restriction in the initial values 

P(0)=L(0), Q(0) = £>(0), (2.16) 
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whose number 2N is the same as that of the Calogero model. The solution of the initial value 
problem of the above equations of motion is given by 

Q(t) = D(0)cosut + ^-Uincjt, Pit) = L(0) cosoot - uD(0) smut. 

v (2.17) 

& X ± (t) = e ±ia '* J L ± (0) 

Substitution of the above solution ()2.17|) into the definition of P(t) and Q(t) in eq. ()2.13|) re- 
produces the solution of the initial value problem of the Calogero model ()2.12j) . In particular, 
the coordinates of the Calogero model is given by the eigenvalues of the matrix-valued harmonic 
oscillator Q(t) in eq. ()2.17|) . That is the essence of the projection method. 

The restrictions on the initial values ()2.16|) can be explained in terms of the constraints on the 
variables Q(t) and P(t). From eqs. ()2.6|) and ([2. 13)1 . one can derive 

[Q(t), P{t)) = U(0, t) [D(t),L(t)]tf(0, t) 
= iaU(0,t)(T - E)rf(0,t). 

By use of the property of the time-evolution unitary matrix (|2.11|) . one obtains 

[Q,P] = ia(T-E), (2.18) 

which poses N(N — 1) constraints on the "unconstrained" variables Q(t) and P(t) given by two 
Hermitian matrices whose off-diagonal elements are pure imaginaries, i.e., 

Qijit) := qui^Sij + iqijit), P^t) := putySy + ipij(t), 



?«(*), Qij{t), Pii{t), Pij{t) E R, q i:j {t) = -qji(t), Pij(t) = -Pji{t). 



The above restriction (J2.19|) is consistent with the solution of the initial value problem (J2.17j) . 
The number of the independent variables in Q(t) and P(t) given by eq. (|2.19J) is A^(A^ + 1). One 
thus reproduces the number of the initial values 2iV as the degrees of freedom of the constrained 
system, iV(iV + 1) — iV(iV — 1) = 2N. The transformation (|2.13|) thus should be interpreted as a 
map of the Calogero model into the harmonic oscillator (|2.15|) of the Hermitian matrix given by 
eq. (J2~T?H) with the constraints (J2~T%|) . 



Ill Integrable Discretization 

As we have discussed in the previous section, the equations of motion (the Lax equation) of 
the Calogero model can be mapped to those of the matrix-valued harmonic oscillator. We thus 
begin with the integrable discretization of the matrix-valued harmonic oscillator, whose difference 
equations of motion are given as 

A+Q n = — {P n +i + Pn)i 

1 2 (3,1) 
A+P n = ~2 LJ (Qn+l + Qn)i 

where Q n and P n are Hermitian matrices whose initial values are fixed as Q$ = D(0) and Pq = L(0) 
so as to relate them with the difference analogue of the Calogero model. The difference equations 
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given above have the same form as those for one- dimensional harmonic oscillator discretized by the 
energy conservation scheme jU] , which is nothing but the implicit midpoint rule giving a symplectic 
integration method of order 2 (see Theorem VI. 3. 4 in ref. PQ). 
As in eq. ()2.13|) . we introduce the new variables 

X± := P n ± koQ n . 

This brings the difference equations of motion into 



A + X± = ±-iu;(X± +1 + X„ d 



which is equivalent to 

(l =F huAt)x± +1 = (l ± hu/\t)xt (3.2) 

Defining the rescaled time-step by 

2 uAt . . 

At := — arctan , (3-3) 

uj 2 

the recursion relation (J3.2|) is rewritten as 

1 _|_ iujAt 

y± 1 - 1 - 2 v± _ iiwAr y± 

+1 1 rn ^At " n ' 

1 + 2 

Thus the solution of the initial value problem of the discrete harmonic oscillator is 



or 



Q n = D(0) cosnuAr H — -^(0) smnuAr, 

uj (3.4) 

P n = L(0) cos no; At — ujD(0) sin no; At, 

in terms of Q n and P n . The above solution ()3.4j) is exactly the same as that for the harmonic 
oscillator with the continuous time (j2.17|) up to time-rescale ()3.3|) . 

The recursion relation ()3.2|) provides us with an efficient way to construct the constants of 
motion of the discrete harmonic oscillator. Consider an arbitrary product of X^~ +1 and X~ +1 like 
(Xn+iY 1 {X~ +l ) mi • • • , for li, mi = 0, 1, 2 • • • , and reverse the time for one discrete time-step using 
the recursion relation ()3.2)1 . Then one obtains 

{X^ +1 ) h {X- +1 ) mi ■■■ = e^X^'e-^^X-r • • • , . 

= e i ^ 1+ -- (mi+ - ))Ar (X+) /l (X7r i 

Taking the trace of the above relation, one obtains 

Tr(IIfti) lj K + i)1 = exp(i^Ar ^{h ~ m k ))Tr(Jl{X+) l '{X;) m >). 
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Thus one can construct constants of motion of the discrete harmonic oscillator in a way parallel 
to what is done for the Calogero model (|2.4j) by 



m i:j 



U Tl (U( X n) k3 ( X n) m -), when J^kj = J2 

i j i,j i,j 

as well as the matrix-valued constants of motion by 

(x+ +1 )'Hx- +1 r • • • = W) mi • • • , 

when l\ + • • • = nil + • ■ ■ ■ In particular, the following quantity, 

[X+, X~\ = -2ua{T - E) [= [X+, X"]] (3.6) 

is conserved since it is a special case of the above matrix-valued constants of motion. Note that 
the constants of motion of the matrix-valued harmonic oscillator with continuous time can be 
given by the same formulas. 

As we have confirmed, the solution of the harmonic oscillator with discrete time (|3.4|) and that 
with the continuous time (|2.17|) agree up to time-rescale (|3.3jl . Thus the eigenvalue of Q n must 
trace the same trajectory of the Calogero model with continuous time. We shall discuss it more 
in detail. 

Since the relations Q n = Q{nAr) and P n = P{nAr) hold as a consequence of eqs. (|2.17|) . (|3.3|) 
and (|3.4j) . we also have analogous relations with ()2.13|) for Q n and P n , 

D n = UlQ n U n = D(nAr), L n = U\P n U n = L(nAr), 

where the unitary matrix U n is also given by the corresponding matrix in the theory for the model 
with continuous time: 

U n = U(0,nAr). (3.8) 
Introduce in analogy with iJ 1 in eq. (|2.3p . i.e., 

L± := L n ± iuD n = UlX±U n (3.9) 
and substitute it into the recursion relation ()3.2|) . Then one obtains 

(l =f l -iuAt)U n+l Lt +l U ] n+1 = (l ± ^xuAt)U n L±Ul 

By multiplying with U\ on the left and U n+ i on the right, one obtains a recursion relation of the 
matrices 



(l =f huAt)S n L± +1 = (l ± kuAty±S n , 



(3.10) 
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which we shall call the discrete Lax equations, or in short, the dLax equations. They will play 
an essential role in the construction of the constants of motion for the discrete time model. From 
the definition of the unitary matrix S n and using the semigroup property of the time evolution 
unitary matrix U(t',t) (|2.10j) . one obtains 

S n := U ] n U n+l = U(nAr, (n+ I) At) = S n [x n , p n , a, u; At), 
U n = S1S2 ■ ■ ■ S n -i, Uq = E. 

This indicates that S n is the one-step time-evolution matrix. We should note that the discrete 
inhomogeneous Lax's equation introduced for the Calogero-Moser model (the Hamiltonian (J1.1J) 
with u = 0) in ref. j2I] inspired us with the above recursion equation (|3.10|) . However, the explicit 
forms of the Lax pair and the derivation of the recursion relation of this paper are different from 
those in ref. |21| . 

The recursion relation ()3.10|) can be interpreted as the discrete time analogue of the Lax 
equation of the Calogero model because of the following two reasons. The first reason is that the 
dLax equations reduces to the Lax equations of the Calogero model (|2.2jl in the continuous 

time limit, At — > 0. The other reason is that the dLax equations also conserve the constants of 
motion of the Calogero model with continuous time. We shall discuss them more in detail. 

From the Taylor expansions of L n+ \ = L((n + l)Ar) and S n = U(t f = nAr,t = (n + l)Ar) 
together with the expression in the formal Dyson series ()2.9|) at t = nAr, one obtains 

AT± 

L- +1 ~L- + ^At + 0(An (312) 
S n ~ E + M n At + 0(At 2 ), M n := M(nAr). 

Substitution of the above expressions into the dLax equations (|3.10|) yields 

(1 t l^At)(E + M n At)(L± + ^At) = (1 ± UuAt)L±(E + M n At). 

Dividing the above relation by At and taking the limit At — > 0, one gets 
H7~± 

<lL " [L±M n ]±koli, 



dt 

which is nothing but the Lax equation of the Calogero model with continuous time. 

In a way parallel to the construction of the constants of motion of the discrete harmonic 
oscillator, one can construct the constants of motion of the dLax equations ()3.10jl . Using the 
definition (j3.9}l of in eq. ()3.5j) and multiplying by U\ and U n+ \ respectively from the left and 
the right, one gets 

S n (L+ +1 )^n+ir • ■ • = e^(Lt) h e- m ^(L-r --'S n 

_ giw(/iH (mi+— ))Ar^+yi ^Tj-^ m i ... 5* 

Thus one has 

(Lt + i) h (L- +i r ■■■ = e-^ + --^ + -» A ^t(L+)^(^) mi ■ ■ ■ S n . 
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Since the trace of an arbitrary product of matrices is invariant under any cyclic change of the 
order of the matrices, one has 

Tr(n^n + i)^+in) = exp(icAr £(Z fc - m k ))Tr (]J(L+) l > '(!£)"*) , 

3 k 3 

which leads to 

i j k,l i j 

One thus concludes that 

X\^{X\{L + n ) l -{L-) m ^) (3.13) 

» i 

gives the constant of motion of the dLax equations (|3.10|) when 

This is in complete agreement with the situation for the Calogero model with continuous time (|2.4J) . 

Since should have the same form as L ± for the continuous time model, the commutator 
between and L~ also should be a constant matrix as in eq. (|2.fij) . In the continuous time 
theory, this constant matrix is associated with the nontrivial constraints of the matrix-valued 
harmonic oscillator ()2.18|) . One can observe the same situations in our discrete time model ()3.10|) . 
Substitution of eq. ()3.9|) into [L+,L~] with the help of eq. ()3.6j) yields 

[L+,L-] = Ul[X+,X-]U n 

= Ul(-2ua(T - E))U n . 

Since U n and also S n are made from the time-evolution unitary matrix U (f, t) as in eqs. ()3.8|) and 
1)3.11]) . they also have a constant eigenvector *[!,...,!] whose eigenvalues are unity, 



U n % . . . , 1] = . . . , 1], [1, . . . , l]U n = [1, . . . , 1] <S> U n T = TU n = T, 
S n % ...,!}= % . . . , 1], [1, . . . , l]S n = [1, . . . , 1] & S n T = TS n = T. 



(3.14) 



Thus one immediately obtains 

[L+,L-] = -2ua(T-E). (3.15) 
using the property ()3.14j) of U n . The relation ()3.15|) is also equivalent to 

[D n ,L n ]=ia(T-E). 

Note that eq. ()3.7j) provides the most general form of L n that satisfies the above relation together 
with the diagonal matrix (D n )ij = x^iy 
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IV Difference Equations of Motion 

Though we have developed an integrable discretization for the Calogero model in the previous 
section, the difference equations of motion we have obtained is nevertheless a formal expression. 
An explicit expression for S n is necessary in order to obtain the dLax equations for the Calogero 
model explicitly. To do this, we restrict the number of particles to two in the discussion of the 
explicit form of the S n matrix. 

From the dLax equations (JXTUJ), one obtains 



J-'n+l 



L n+1 



1 ± hcoAt . , 

2 Ot T ± O 

1 =F puAt 



Solution of the above equations with respect to D n+1 and L n+ i are given by 
1 



D 



n+l 



_Qt 
2iuT" 



1 + + 1 



I - \iuAt 



\iuoAt 
-L~ 



S n , 



2 n 



l + >At 
1 McuAt 



1 



±ia;At 
-K 



2 1 2 

Since D n+ i is diagonal, diagonalizes 

1 



l + \\ujAt n 



(4.1) 



1 + iiwAi r 



1'iujAt' 



kuAt 

LI 



1 + \iuAt 



S n thus can be constructed from the eigenvectors v^ n , % — 1, • • • , N of the above matrix. Normal- 
ized eigenvectors s^ n that constitute = [s 1<n ■ ■ ■ s N)n ] should have the following form, 



V; 



N (v- )■ 



N 

B 



(4.2) 



so as to satisfy the condition ([3. 14)1 . Ordering of the eigenvectors s i n is uniquely determined by 
the Taylor expansion of S n (|3.12j) . or equivalently 



lim S n = E. 

Considering the eqs. (J4.1j) (|4.2|) and (J4.3j) all together, one confirms that 

2aiAt 



(4.3) 



S n 



E 



Y n - 2aiAt + ^Aa 2 At 2 + F n 2 

1 - 



1 -1 
-1 1 



E + {\M- hn At - M r>n Ar 



1 



(4.4) 



E + M n At (At -> 0) 



gives the explicit form of the 5 n matrix of the dLax equations for the Calogero model for the 
two-body case. It is straightforward to verify that the dLax equations ()3.10|) with the above S n 
matrix f)4.4j) yields the difference equations of motion ()1.2|) shown in the introduction. 

In order to confirm how well the dLax equations (|3.1U|) with the S n matrix (J4.4j) or the difference 
equations of motion (|1.2j) of the Calogero model with discrete time trace the behavior of that with 
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continuous time, the explicit expression of the analytic solution for the continuous-time model is of 
great help, which is given by the eigenvalues of Q(t) in eq. ([2.17)1 . But there is another derivation. 
Since U(t',t) is related to S n according to eq. ([3.1 1|) with t' := nAr and t := (n + l)Ar, one can 
obtain the explicit form of U(t',t) for the two-body case from that of S n ()4.4[) : 



U(t',t) = E + 



2ai sin oj(t-t') 



y 



2ai sinoj(t— t') 



4a 2 sin 2 u(t-t') 



+ y 2 



i -i 
i i 



(4.5) 



y : = (xi - x 2 ) 2 cosuj(t - t') + (pi — p 2 )(xi - x 2 ) 



sin uj(t — if) 



UJ 



Substitution of the explicit forms of U(0, t) (|4.5|) and -f/ ± (0) for the two-body case into eq. (|2.12j) . 
one obtains the explicit form of the analytic solution of the two-body Calogero model with con- 
tinuous time: 



xiit) 



(xin + x 2 ,o) cos ujt H — (pi o + p 2 ,o) sinut 

UJ 

1 / 4a 2 sin 2 ujt 



2(xi j0 - x 2fi 
=:xi[a?o,P ; a,uj;t] 



+ 



UJ" 



(x lfi - x 2fi ) 2 cos ujt + (px o - p 2 ,o)(xi,o - x 2 ,o) 



sin utl 2 



CO 



=2 (Pi,o+P2,o) coswt - ^(^1,0 + x 2 ,o) sinut 



+ 



1 



2(* 1 ,o-x 2 ,oh/^^ + 



2a 2 sin 2ujt 



CO 



Ol,0 - X2,o) 2 COS Cttf + (pifi - p 2 ,o) (^1,0 - X2,o) 

sin cut" 



(4.6) 



(xi )0 - x 2 ,o) cos u;i + (pi >0 - p 2 ,o) (^1,0 - x 2)0 ) - 



X 



-w(a;i o - £2,o) 2 sin cut + (p 1)0 - p 2 ,o) (^i,o - x 2)0 ) cos cut 



^[sbcPo? ^;*] 



x 2 (t) =x 1 (t) 



1^2 



P2(t) = Pl(t) 



l<->2 



The above expressions are used to numerically display the exact analytic results in the following 
figures. The initial values, the interaction parameter, the strength of the external harmonic well 
are set at (xi,pi,x 2 ,p 2 )\ t=0 = (—4.00,5.00,2.00,1.00), a = 3.00 and uj = 0.314 throughout the 
numerical calculation in the following. 

Figure Q presents the time evolution of the coordinates and the momenta generated by the 
analytic solution and the difference equations of motion (jl.2J) that gives the super-integrable 
discretization. The time-step in the super-integrable discretization is set at At = 1.00. The 
relation between the time t and the number of iterations n is given by n = t/ At where Ar = 
0.991903. The time interval 240 x ^ = 4802 < t < 4842 = 242 x ^ corresponds to 4841 < n < 
4881 in terms of the number of iterations for the discrete case. As has been confirmed by the 
correspondence between the matrices L and D for the super-integrable discretization and those 
for the continuous time model ()3.7|h the solution of the discrete equations of motion ()1.2|) gives a 
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sequence of canonical variables that are "sampled" from the orbit of the exact analytic solution, 
i.e., 

Xi, n = Xi(t = nAr), p i)U = pi(t = nAr). (4.7) 
Even though the coordinates and the momenta generated by the difference equations of mo- 



8 15 




t 4 

CD 

a 

o 

a 2 



4810 4820 4830 4840 

time t 









f- 































4810 4820 4830 
particle 1 particle 2 discrete 



4840 
time t 
analytic 



Figure 1: The time evolution of the coordinates and the momenta generated by the analytic 
solution and the difference equations of motion (jl.2|) . 

tion (jl.2j) are computed through sufficiently large numbers of iterations, they are in quite good 
agreement with those of the exact analytic solution, which gives a good numerical confirmation 
of the correspondence (|4.7j) . 

One can confirm the precise agreement of the two solutions in a more reliable manner by 
observing the relative errors of the constants of motion. Using eq. (|2.4J) (or eq. (|3.13j) ). we have 
already given three constants of motion of the Calogero model (|2.5j) . Note that C\ and I\ cor- 
respond to the Hamiltonian for the center of mass and the Calogero Hamiltonian, respectively. 
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Figure 2: The relative errors of the constants of motion Cf := Cj^/C^o, n = t/Ar, of the discrete 
solutions given by the super- integrable discretization (|1.2|) . 

Those three constants of motion give a set of independent conserved quantities. But just for 
convenience, we introduce a set of slightly modified constants of motion: 

Ci = (pi+p 2 ) 2 + uj 2 (x 1 +x 2 )\ 

4a 2 

C 2 := 2/i - Ci = (pi - p 2 ) 2 + ^ 2 (a;i - £2) 2 + 7 r», /„ s x 

(xi - x 2 y (.4.8J 

C 2 - J 2 2a 2 (x 2 + a; 2 ) 

<^3 := , 2 = O1P2 - X2P1) + "7 xo - - 

4a;^ (xi — X2) 

Note that C2 is the Hamiltonian for the relative coordinates and C3 can be interpreted as a "mod- 
ified quadratic angular momentum." We use the above Cj's, % — 1, 2, 3, for numerical calculation. 

Figure El presents the relative errors of the constants of motion Cf := Ci,„/Cj,o, n = t/Ar of 
the discrete solutions given by the super-integrable discretization (jl.2J) . The offsets (= (i— 1) x 1.0) 
are added for convenience of presentation. Even though one sees the growth of errors caused by 
round-off errors that is inevitable in any numerical analysis, one confirms that all the constants of 
motion are conserved with a precision of the machine epsilon times the number of iterations. Thus 
one concludes that fig. El presents a direct numerical verification of the fact that the difference 
equation of motion ()1.2j) preserves the super-integrability of the Calogero model. 

Let us compare the above results with numerical demonstrations by two other discretizations. 
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The first one is the energy conservation scheme (see pp. 159-161 in ref. pQ), 

a ^2/ \ V - ^ 0, {{Xi^+Y Xj^ n +i) -\~ {Xi^ n Xj^^j (4-9) 
^■+Pi,n = —-^> {%i,n+l + %i,n) + 7 ' ~ \2? ~ ^2 ' 

which keeps two constants of motion C\ and l\ (consequently, C2) exactly for arbitrary number 
of particles N. As one can see, the scheme (14. 9 [I is an implicit scheme, which involves numerical 
solution of simultaneous algebraic equations. The other scheme is the symplectic Euler method (see 
Theorem VI. 3. 3 in ref. pQ), 

Ajf-Xi^n Pi,nj 

N 2 

A +Phn = -u 2 x i>n+1 + 7 T W ^ 4 ' 10 ^ 

which is an explicit scheme. Note that A + denotes the advanced time-difference (J1.3)) . 

Figures El and |U present the relative errors C^ c ' se := Cj >n /Cj j0 , n = t/At of the discrete solution 
given by the energy conservation scheme (ec) and the symplectic Euler method (se) for the two 
body case (N = 2). We should note that we do not present all the data of each relative errors 
to keep the size of the data file reasonable and that the offsets (= (i — 1) x 1.0) are added for 
convenience of presentation. The initial condition, the coupling parameter and the strength of the 
harmonic confinement are the same as those given for the numerical calculation that gives figs. Q 
and El The discrete time step is set at At = 0.200 for these schemes. 

One observes that the relative errors of the energies (Cf c and C| c ) calculated by the energy 
conservation scheme ()4.9|) remain within the order of the machine epsilon (= 10~ 16 ) times the 
number of iteration (= t/At), which gives a numerical confirmation that the energies are conserved 
by the scheme. On the other hand, one also confirms that the relative errors (except for C\ c and 
C^) are much larger than the order of machine epsilon times the number of iteration, even though 
one cannot observe growth of the relative errors in both two figures. These errors are not brought 
about by the round-off errors at the order of machine epsilon times the number of iteration, but 
by the schemes themselves. One perceives that there is a clear distinction between fig. |21 or the 
super-integrable discretization (jl.2j) and figs. El and 0] or the energy conservation scheme (J4.9|) and 
the symplectic Euler method (J4.10)) . 

One of the characteristics of the super-integrable system is that its bounded orbit is always 
closed in the phase space and, in particular, in the configuration space. As presented in fig. 
the orbit in the xi-x% plane, or the "Lissajous plot" in other words, thus provides the distinctest 
way of comparing the three discretizations. The time interval of the Lissajous plots is < t < 
200, which corresponds to ten periods of motion. While the orbits generated by the energy 
conservation scheme and the symplectic Euler method are not closed, the orbit by the super- 
integrable discretization is always on that of the exact analytic solution, which is, of course, 
closed. 
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Figure 3: The relative errors of the constants of motion, Cf c := C^/C^o, n = t/At, of the discrete 
solution given by the energy conservation scheme (J4.9j) . 
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Figure 4: The relative errors of the constants of motion, Cf e := Cj^/C^o, n = t/ At, of the discrete 
solution given by the symplectic Euler method ()4.1()|) . 
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Figure 5: The "Lissajous plots" 



V Summary and Concluding Remarks 
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The aim of the paper was to present a time-discretization for the Calogero model that maintains its 
super-integrable structure. As discussed in section [TT1 the super-integrable structure of the original 
Calogero model is provided by the Lax formulation and the projection method, which are closely 
related to the equations of motion of the harmonic oscillator. With the help of the integrable 
discretization of the harmonic oscillator jjjj, the Lax formulation and the projection method are 
discretized. As a consequence, a time discretization that preserves the super-integrability of the 
Calogero model is presented in section ITTT1 In particular for the two-body case, an explicit form of 
the difference equations of motion (jl.2j) is obtained in section HVl The difference equations give an 
explicit scheme for time-integration. Numerical results by the super-integrable discretization (jl.2J) 
together with comparison with those by the energy conservation scheme ()4.9j) as well as the 
symplectic Euler method (j4.10)) are presented in five figures in section IIV| which give a intuitive 
numerical verification of the super-integrability of the scheme. 

Lastly, we should give several remarks on previous studies that are relevant to the present 
work. In refs. j2U [221 12H] ( an d also in ref. 0), integrable discretizations of the rational Calogero- 
Moser model (the case u = 0) and its trigonometric, elliptic and "relativistic" (a g-difference 
generalization with respect to space coordinates) generalizations were presented. The structure of 
the projection method also underlies these integrable discretizations, but the interaction parameter 
of the continuous-time models and the time-step of the discrete-time models are related with each 
other in the integrable discretizations above. On the other hand, our discretization preserves not 
only the integrablity but also the super-integrablity of the Calogero model and both the time-step 
At and the interaction parameter a independently appear in the discrete equations of motion. 
Thus our discretization is apparently different from that of the previous studies. The possibility 
of the mutual penetration of the particles in the time-discrete model was reported in ref. [21] , which 
is certainly unlike the continuous-time case. But the scheme given there is implicit and it could 
not privide a way to verify this possiblity. In our super-integrable discretization of the present 
work, however, we gave an explicit scheme of the Calogero model for the two-body case (jl.2|) and 
numerically observed in fig. ^ that there was no penetration of the particles. And this observation 
should be the same for the Calogero-Moser case corresponding to the limit uj — > of the present 
work. The comparison of the two different discretizations together with additional consideration 
on the super-integrability of the rational Calogero-Moser model [21] will be presented in a separate 
paper. Further studies on the trigonometric, elliptic and relativistic Calogero-Moser model along 
the line of our super-integrable discretization of the Calogero model as well as to obtain explicit 
forms of the difference equations of motion for general TV-body case (or at least 3-body case) is 
worthy of interest. 
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